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We investigate the long-standing problem of the hole localization at the A1 impurity in quartz Si02, using a 
relatively recent DFT hybrid-functional method in which the exchange fraction is obtained ab initio, based on 
an analogy with the static many-body COHSEX approximation to the electron self-energy. As the amount of 
the admixed exact exchange in hybrid functionals has been shown to be determinant for properly capturing 
the hole localization, this problem constitutes a prototypical benchmark for the accuracy of the method, 
allowing one to assess to what extent self-interaction effects are avoided. We obtain good results in terms of 
description of the charge localization and structural distortion around the A1 center, improving with respect 
to the more popular B3LYP hybrid-functional approach. We also discuss the accuracy of computed hyperfine 
parameters, by comparison with previous calculations based on other self-interaction-free methods, as well as 
experimental values. We discuss and rationalize the limitations of our approach in computing defect-related 
excitation energies in low-dielectric-constant insulators. 


I. INTRODUCTION 

The A1 impurity is one of the most commonly observed 
defects in irradiated quartz Si 02 , by which a tetrava- 
lent Si cation is replaced with a trivalent A1 atom. The 
Al/Si substitution results in an unpaired electron (or, 
in an equivalent description, a hole) which, based on 
early experimental observations, has been identified to 
be trapped in a nonbonding 2p orbital of an O atom 
surrounding the substitutional Alfii^ The corresponding 
neutral A1 defect center (hereafter also denoted [AlO^]^) 
is magnetically active, and has been the subject of ex¬ 
tensive characterization by electron paramagnetic res¬ 
onance (EPR)j^“— also in combination with absorption 
spectroscopy i^“— 

Prom the theory side, it has been recognized that the 
reproduction of the experimentally observed features of 
the [A 104 ]°center is subordinate to a correct description 
of the hole localization properties. Since early density- 
functional theory (DFT) calculations using (semi)local 
density functionals gave a wrong picture, predicting the 
hole to be delocalized over all the four A1 nearest- 
neighbor O atomS)^^— this problem has been identified as 
a challenging testing ground for novel density-functional 
methods 

The wrong description of (semi)local DFT functionals 
has been ascribed to incomplete cancellation of the self¬ 
interaction (SI) brought in by the Hartree termjiid^ since 
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calculations using exactly Sl-free Hamiltonians, such as 
unrestricted Hartree-Fock (UHF) and self-interaction- 
corrected (SIC) DFT, yield the correct hole localiza¬ 
tion. Pacchioni et al. also performed unrestricted 
second-order Mpller-Plesset perturbation theory (UMP2) 
calculations?^ concluding that no appreciable role is 
played by correlation in the specific problem at hand. 

A now routinely used approach allowing to par¬ 
tially amend the SI error is based on admixing a frac¬ 
tion of Hartree-Fock exact exchange (EXX) to semilo¬ 
cal exchange-correlation functionals. However, popular 
recipes for such hybrid functionals, prescribing 20% or 
25% of EXX, such as in B3LYP,iid^ PBE0?i^ and HSE06 
(Ref. [ 13 ) have proven to fail in describing charge localiza¬ 
tion in Al-doped silicaj ^^d^d^~ — Thus, it has been argued 
that a large amount of EXX would be necessary to ob¬ 
tain agreement with experiments much like a large 
enough Hubbard U parameter is needed within DFT-I-H 
to solve the same problem,™ Eor instance. To et al. 
found that a semi-empirical hybrid functional including 
42% of EXX (called BBIK functional) yields the correct 
picture)^ 

In view of the preceding work, it may be concluded that 
some amount of empiricism is required to tackle the Al- 
impurity problem within DET, casting doubts on its ac¬ 
tual predictive power. Recently, a rationale for the value 
that the EXX fraction takes in extended systems has 
been put forwardi^^i^i Starting from the many-body GW 
approximation to the electron self-energy E(a;), one ob¬ 
tains, in its static limit (w —)• 0), the so-called Coulomb- 
hole-plus-screened-exchange (COHSEX) approximation 
to By performing a spatial average of the electron 
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gas polarization function, one is finally left with a simple 
analytical expression for the EXX fraction, which is now 
expressed in terms of the macroscopic electronic dielectric 
constant of the material. Thus, the EXX fraction can be 
obtained ab initio for a given material once its dielectric 
constant is computed within DFT; a self-consistent ap¬ 
proach to the definition of such self-consistent dielectric- 
dependent hybrid functional has been recently proposed 
and tested for various properties of oxide semiconduc¬ 
tors and insulators Notice that COHSEX treats the 
exchange term exactly, and hence is SI free; it thus con¬ 
stitutes the ideal candidate for modeling systems whose 
ground state is dominated by the classical Hartree and 
electron exchange Fock interactions. 


II. COMPUTATIONAL DETAILS 

All-electron DFT calculations were performed within 
the linear combination of atomic orbitals approach as 
implemented in the CRYSTAl09 code i^^’^^ The follow¬ 
ing all-electron basis sets were em ploy ed: 66-21G* for 
Si (Ref. 113 ), 8-411(dl) for O (Ref.lj, 85-llG* for A1 
(Ref.m. A full-range hybrid functional^ was adopted 
for the treatment of exchange and correlation; the frac¬ 
tion of admixed EXX (a) was evaluated based on the 
above-mentioned relationship with the macroscopic elec¬ 
tronic dielectric constant (coo), ct ~ l/coo,— and ob¬ 
tained self-consistently for pristine quartz Si 02 , following 
the procedure illustrated in Ref. [23- The resulting func¬ 
tional is referred to as “sc-PBEOoe,,^” in the following, as 
it is de facto a self-consistent (sc) re-parametrization of 
the PBEO hybrid functionali^^ The dielectric constant Coo 
was computed within the coupled-perturbed Kohn-Sham 
method implemented in the CRYSTAlOQ code.— 

The A1 center was modeled in an embedding 2x2x2 
quartz Si 02 supercell (72 atoms) with the atomic posi¬ 
tions and lattice parameters fully relaxed for the bulk cell 
using the sc-PBEOoe,,^ functional. For the defective su¬ 
percell, further optimization of the atomic positions was 
carried out at hxed lattice parameters. The standard 
thresholds defined in CRYSTAl09 were adopted in all ge¬ 
ometry optimizations.— The Brillouin zone was sampled 
by using 8 fc-points in the irreducible wedge— 

The charge-transition levels were computed according 
to the formalism illustrated in Refs. [13 andjl^- In partic¬ 
ular, total-energy differences relative to the defect charge 
state variation were computed using defect Kohn-Sham 
(KS) eigenvalues, following the approach proposed in 
Ref. 113. The Is KS eigenvalue of Si was taken as ref¬ 
erence for aligning band structures in defect and bulk 
calculations. The spurious electrostatic interaction be¬ 
tween image charged defects was accounted for by cor¬ 
recting the KS eigenvalues according to the procedure 
illustrated by Chen and Pasquarello^ and based on the 
Makov-Payne correction scheme— 


III. RESULTS AND DISCUSSION 

Within the sc-PBEOoe,,^ approach, the dielectric con¬ 
stant of quartz Si 02 turns out to be Coo = 2.15 (the 
experimental value is 2.38, see Ref. [13 and references 
therein), which corresponds to an exchange fraction a = 
46.5%. The method predicts a somewhat overestimated 
band gap of 11.6 eV, whereas various experiments mea¬ 
sured it in the quite broad range of 8 — 10 eV— The 
failure of the sc-PBEOa,,,,^ functional in computing reli¬ 
able band gaps for some insulators with very low dielec¬ 
tric constants (and correspondingly large band gaps) was 
already reported in Refs. 113 and 113. being particularly 
serious when the geometry is re-optimized at each self- 
consistency step (see for instance the case of MgO, which 
has a dielectric constant of ~ 3) — 

However, here we are mainly concerned with the de¬ 
scription of the hole localization at the [AlO^J^center, 
which is a ground state property; the related spectro¬ 
scopic features will be deferred to a separate discussion 
in Section IIII Cl 


A. Structural deformation and hole localization 

Figure [IKa) shows the local atomic structure of Si 02 
around the A1 impurity, as found by minimizing the 
total electronic energy using the sc-PBEOa,,,,^ functional 
and allowing for symmetry-breaking atomic relaxations, 
which amounts at independently optimizing the positions 
of each atom of the (Oi, O 2 ) and (O 3 , O 4 ) oxygen pairs 
(where the two O within each pair are equivalent to each 
other in the bulk Si 02 structure). The emerging pic¬ 
ture agrees very well with that obtained within other 
rigorously Sl-free approaches, such as UHF— and SIC- 
DFT— In particular, the hole introduced by the substi¬ 
tutional A1 atom is trapped at the 0(1) atom, its wave- 
function exhibiting purely 2p character, with the corre¬ 
sponding orbital lying almost perpendicularly to the Al- 
0-Si plane. As a consequence of the charge localization, 
the local atomic structure distorts considerably: the 0 ( 1 ) 
atom moves away from the A1 atom, resulting in an av¬ 
erage equilibrium Al-0 distance 13% larger compared to 
the other Al-0 distances (see also Table IJ). 

The results of our calculations are in agreement 
with the experimental evidence as obtained from EPR 
investigations^ (i) the hole localizes in the nonbonding 
2p orbital of the O atom corresponding to the longer Si- 
0 -type bond in pure Si 02 ; (ii) the 2p orbital is perpen¬ 
dicular to the Al-O-Si plane; (hi) the localization causes 
the hole-bearing oxygen to move away 12 % farther from 
the A1 center with respect to the other O atoms. 

For comparison, we also tested the performance of the 
Becke’s three-parameter hybrid functional (B3LYP)^d^ 
incorporating 20% of EXX. From Table |T] it is inferred 
that, while some asymmetry is still present in the result¬ 
ing optimized structure, the elongation of the Al-O(l) 
distance is now at most 6 % larger than the other Al-0 
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TABLE I. Nearest-neighbor Si-O and Al-O distances (A) for 
pure bulk and Al-doped Si02 ([A 104 ]°), respectively. The O 
atoms are labeled according to Figure [T] 


Functional 

Pure Si02 
sc-PBEOqc^^ 


[A104]° 


SC-PBEOOfe,, 

, B3LYP^ 

B3LYP^ 

0(1) 

1.620 

1.910 

1.809 

1.826 

0(2) 

1.620 

1.699 

1.758 

1.749 

0(3) 

1.616 

1.687 

1.705 

1.705 

0(4) 

1.616 

1.691 

1.700 

1.700 


Starting geometry for optimization: pure Si02 structure. 
Starting geometry for optimization: sc-PBEOoej,^-optimized 
[A 104 ]'^ structure. 


(a) (b) 

FIG. 1. Ball-and-stick representation (Si, O and A1 atoms 
are shown as yellow, red, and pink spheres, respectively) of 
the local atomic structure around the A1 impurity, as ob¬ 
tained from geometry optimization at the (a) sc-PBE0Qe„„, 
and (b) B3LYP level of theory. Isosurface of the spin den¬ 
sity associated with the unpaired electron introduced by the 
[A 104 ]°center is shown. For bond distances see Table H] 

distances. The Mulliken population analysis presented 
in Table m shows that in the B3LYP ground state the 
unpaired electron charge density is distributed over the 
(0(1), 0(2)) pair, with a substantial contribution also 
from a 2p orbital of the 0(2) atom [see also Figure HKb)], 
at variance with the experimental evidence. The final 
picture is not dissimilar if the optimization is started 
from the distorted structure of the [A 104 ]°center opti¬ 
mized at the sc-PBEOag^leveli^ the Al-O(l) elongation 
now amounts at about 7%, and the unpaired electron is 
delocalized over the (0(1), 0(2)) pair with roughly the 
same proportions as reported in Table HU In some of the 
previous B3LYP studies of Al-doped Si02 the inequiva¬ 
lence of the 0(1) and 0(2) sites is not even qualitatively 
captured (same Al-O distances), and the hole is found 
evenly delocalized over two, or even over all the four O 
atoms nearest to the A1 impurityj ^^d^ 


B. EPR parameters 

In order to further confirm the better performance 
of the sc-PBEOae^ approach with respect to B3LYP, 


TABLE 11. Spin population of the O atoms belonging to the 
[A 104 ]°center, and EPR hyperfine parameters of the hole¬ 
bearing and of the ^^Al. The optimized structures are 
obtained starting from ideal bulk Si 02 . 


Functional 

Geometry 

SC-PBEOOe,, 

SC-PBEOOe,, 

„ B3LYP 
„ B3LYP 

B3LYP 

SC-PBEOCTe,,^ 

Expt. 

(Refs. ^ and^ 



Spin population 


0(1) 

0.95 

0.58 

0.81 


0(2) 

0.02 

0.28 

0.07 


0(3) 

< 0.01 

0.01 

0.03 


0(4) 

< 0.01 

0.04 

< 0.01 



'"0(1) 

hyperfine matrix (G)^ 


-^iso 

-42.6 

-26.1 

-34.7 

-26.0 


-94.3 

-61.7 

-83.7 

-85.0 

B 2 

47.1 

30.7 

41.8 

41.2 

B 3 

47.2 

31.0 

41.9 

43.8 


^^Al hyperfine matrix (G)'‘ 



-5.0 

-8.1 

-5.4 

-5.8 

Bi 

-0.4 

-0.2 

-0.4 

-0.4 

B 2 

-0.4 

-0.1 

-0.3 

-0.3 

B 3 

0.8 

0.3 

0.6 

0.7 


^ The principal values of the anisotropic hyperfine matrix are 
listed so that < B 2 < B 3 . 


we computed hyperfine parameters relative to the hole¬ 
bearing 0(1) atom. The hyperfine coupling matrix, de¬ 
scribing the magnetic interaction of the spin of the un¬ 
paired electron with the spin of the neighboring nuclei 
(^^O and ^^Al), is conveniently divided into an isotropic 
(spherically symmetric) and an anisotropic (dipolar) 
part. The isotropic part (denoted Aiso) is proportional to 
the electron spin density at the nucleus, and, as such, the 
dominant contribution to it is caused by spin-polarization 
of the s electrons. Instead, the anisotropic part is related 
to the spin population of orbitals with higher angular 
momentum components; this contribution is commonly 
expressed in terms of a matrix with principal values Bi, 
B 2 and B 3 reported in Table HD Since the hole wave- 
function has purely 2 p character, the dipolar part gives 
direct access to the corresponding spin distribution. In¬ 
stead, the isotropic contribution is notably harder to be 
reproduced, being extremely sensitive to the details of 
the calculation in general, and to the choice of the basis 
set in particulars^ We nevertheless report on it as well 
for the sake of completeness. 

In Table HD the hyperfine parameters are reported for 
the ^^0(1) and ^^Al nuclei. Concerning the anisotropic 
parameters of our sc-PBEOae^ calculations nicely 
capture the experimentally observed strong anisotropy 
along the three axes, and numerical values are in good 
agreement with both experiments and the results of pre¬ 
vious investigations based on Sl-free approaches (UHF, 
UMP2, SIC-DFT)JSS3i46 contrast, B3LYP yields a 
quantitatively wrong picture, the computed parameters 
being substantially smaller in absolute value than exper- 
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TABLE III. Vertical excitation energies (in eV) associated to the [A 104 ]'^center computed with the sc-PBEOoe^ and B3LYP 
functionals at different optimized geometries. Comparison with theoretical results from the literature obtained within time- 
dependent DFT (TDDFT) using the B3LYP and BBIK exchange-correlation approximations, as well as within the outer-valence 
Green’s function (OVGF) approach. Experimental position of the main absorption peak is also reported. 


Functional 

This work 
sc-PBE0ae„„ 

B3LYP 

TDDFT-BBIK 

TDDFT-B3LYP 

OVGF 

Expt. 

Geometry 

sc-PBE0ae„„ 

sc-PBE0ae„„ 

(Ref. 22J 

(Ref. ^ 

(Ref. 20) 

(Ref. Ref.^ 


4.73 

2.91 

3.03 

1.72 

2.74 

2.9, 2.85 


imental ones. The situation is quite the opposite for the 
isotropic part, for which sc-PBE0q;£„^ overestimates the 
absolute value of Aiso, while B3LYP yields it exception¬ 
ally close to experiment. We argue that the latter result 
is fortuitous, in the sense that it is not concomitant with 
a correspondingly more accurate description of the hole 
localization. Firstly, as already discussed, at the B3LYP 
level the 0(2) atom carries a substantial part of the hole- 
related spin density (see Table |TT]) , in disagreement with 
experiment, and accordingly the ^^0(2) EPR parameters 
are of the same order of magnitude as for ^^0(l)f^ in¬ 
stead, when the hole localization is correctly captured, 
such as at the sc-PBEOae„^ level, the former are at least 
an order of magnitude smaller than the latter (see also 
Ref.lH). Secondly, the improvement of the isotropic part 
does not come along with a similar improvement of the 
anisotropic one, which indeed is related to the proper 
description of the 2 p hole wavefunction. 

The above conclusion is also supported by the 
computed superhyperfine matrix of ^^Al: the 

sc-PBEOae„^ functional gives results in quantita¬ 
tive agreement with both experiment and UHF 
calculationswhereas this is not true for B3LYP. 
Notice that the superhyperfine interaction with a dopant 
element like A1 that introduces a hole in the structure is 
usually the only accessible information. In fact, in order 
to measure the O hyperfine constants, enriched 
samples have to be prepared with complex and costly 
procedures 


C. Optical properties of the Al impurity 

The Al impurity in quartz silica has been observed to 
act as a color center, endowing it with the typical smoky 
coloration. However, considerable controversy arose as 
to which absorption feature had to be correlated with 
such observation!^ It was finally concluded that an ab¬ 
sorption peak at about 2.9 eV is to be associated with 
the presence of Al centers and, thus, with the smoky 
coloring!^ 

The observed optical transition should be related to ex¬ 
citation of the hole trapped at the neutral [AlO^J^center 
{q = 0) into the VB, leading to a negatively charged 
defect {q' = —1, [AlO^]”^center); accordingly, from the 
theory side, the optical transition level (0/ — 1) is the 
relevant quantity to be compared with experiments. The 


sc-PBE0ae„^ predicts the computed level to be 4.7 eV 
above the valence band (VB) maximum, nearly 2 eV 
higher in the band gap than in experiments. We attribute 
this disagreement to the already mentioned overestima¬ 
tion of the electronic band gap of bulk quartz provided by 
sc-PBE0ae„^. Consequently, the Al-related defect level is 
wrongly positioned with respect to the VB (^ 5 eV above 
its edge) and this eventually gives rise to the observed 
overestimation of the optical transition energy. 

As a partial workaround, we computed the electronic 
structure using the B3LYP functional, which yields a 
band gap of 8.6 eV for bulk quartz, falling in the range 
of the experimental values; the [A 104 ]°geometry obtained 
within sc-PBEOae„^was instead retained. The spin popu¬ 
lation and EPR parameters computed following this ap¬ 
proach are reported in Table El The hole is again local¬ 
ized on the 0(1) atom, although with a Mulliken density 
lower than that obtained by performing calculations fully 
within the sc-PBE0ae„^scheme. Surprisingly, the com¬ 
puted hyperfine parameters are in even better quantita¬ 
tive agreement with experiments, as the lower spin den¬ 
sity counterbalances the overestimation yielded by the 
sc-PBE0Q!e„^; a similar trend was noticed in the previous 
hybrid-functional investigation of To et al^ The charge 
density distribution analysis for the negatively-charged 
center led us to conclude that the same qualitative pic¬ 
ture is obtained at the sc-PBEOa,;,,^ and B3LYP levels, 
provided that the geometry is kept fixed to the one opti¬ 
mized within sc-PBEOae,,^: the one-particle state corre¬ 
sponding to vertical excitation of the hole to the VB is 
still contributed by the 2p orbitals of the 0(1) atom. The 
computed optical level for such transition is positioned at 
^2.9 eV above the top of the VB, in excellent agreement 
with experiment. For comparison, we report in Table El 
representative results from previous theoretical studies 
for the computed vertical transition energy correlating 
with the experimentally found absorption band with a 
maximum at 2.9 eV. 


IV. CONCLUSIONS 

We have re-investigated the long-standing problem of 
hole localization in Al-doped quartz Si02 using a re¬ 
cently proposed hybrid-functional method in which the 
exchange fraction is consistently determined based on the 
analogy with the many-body COHSEX approximation 
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to the electron self-energy, without empirical prescrip¬ 
tions or fitting to the experimental data. The COH- 
SEX scheme is rigorously Sl-free, and thus constitutes 
the ideal starting point for studying systems in which in¬ 
complete cancellation of SI leads to a qualitatively wrong 
ground state. The neutral A1 impurity in quartz silica is 
just a paradigmatic case: (semi)local or standard hybrid 
DFT functionals fail in capturing the experimentally ev¬ 
idenced hole localization at one of the Al-coordinated O 
atomsi ^^d^ In particular, the failure of popular hybrid 
functionals, such as B3LYP, has been attributed to the 
insufficient amount of EXX admixed. 

The sc-PBEOae,,^ approach allows one to evaluate the 
exchange fraction from first-principles, based on a simple 
relationship with the macroscopic dielectric constant of 
the material. In the case of quartz silica, we obtained it 
to be ~ 46%. The resulting hybrid functional correctly 
reproduces the hole localization at a single O atom sur¬ 
rounding the A1 impurity, also giving an accurate descrip¬ 
tion of the structural distortion around it, and allowing 
to compute EPR parameters in agreement with previ¬ 
ous STfree calculation a^^d^i^^ and experiments^- How¬ 
ever, the defect-related optical spectroscopic features are 
not well-reproduced by the sc-PBEOoe,,^ method. We at¬ 
tribute this failure to the observed overestimation of the 
bulk quartz silica band gap. Using B3LYP on top of the 
sc-PBEOae^ optimized geometry yields a band structure 
in better agreement with experiment and corroborates 
this hypothesis. 

In conclusion, we have shown that the hybrid- 
functional approach tested in this work is capable of 
correcting most of the SI error inherent to (semi)local, 
as well as more popular hybrid, DFT functionals. This 
feature is crucial for adequately describing the ground 
state of defective oxide materials. As far as defect-related 
excitation energies are concerned, their determination 
is subordinate to an accurate calculation of the elec¬ 
tronic structure of the bulk material. This is not always 
the case for low-dielectric-constant insulators in general, 
and for quartz Si02 in particular. However, based on 
our experience, the sc-PBEOofe^ method is able to repro¬ 
duce the whole experimental scenario when point defects 
in moderate gap metal-oxide semiconductors (dielectric 
constants ~ 4 — 6) are addressed. 
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